Halo independent comparison of direct dark matter detection data 
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We extend the halo-independent method of Fox, Liu, and Weiner to include energy resolution and 
efficiency with arbitrary energy dependence, making it more suitable for experiments to use in pre- 
senting their results. Then we compare measurements and upper limits on the direct detection of low 
mass (~10 GeV) weakly interacting massive particles (WIMPs) with spin-independent interactions, 
including the preliminary upper limit on the annual modulation amplitude from the CDMS collab- 
oration. We find that isospin-symmetric couplings are severely constrained, but isospin-violating 
couplings are still possible if for example the local Galactic escape speed is small, as found in recent 
surveys. 
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The nature of dark matter is one of the fundamental 
problems of physics and cosmology. Weakly interacting 
massive particles (WIMPs), i.e. particles with weakly in- 
teracting cross sections and masses in the GeV-10 TeV 
range, are among the best motivated candidates for dark 
matter. Of particular interest is a low mass region, ~ 10 
GeV, suggested by data from three direct dark matter 
experiments: DAMA [I], CoGeNT 013] and CRESST- 
II g]. DAMA and CoGeNT report annual modula- 
tions with the expected characteristics of a WIMP sig- 
nal [S]. CRESST-II observes an excess of events above 
their known background, excess which may be inter- 
preted as due to dark matter WIMPs. Stringent upper 
limits have been placed on dark matter WIMPs by other 
direct detection experiments. The most stringent lim- 
its in the region of ~ 10 GeV WIMPs come from the 
XENONIO 0, XENONIOO [7], SIMPLE [gj, and CDMS 
experiments [5]. All but one of these limits result from 
an upper bound on the total unmodulated event rate. 
The exception is a very recent preliminary result by the 
CDMS collaboration [10], which has searched for an an- 
nual modulation in their data and, not finding it, has 
placed a stringent upper limit on its amplitude. 

In this paper, we compare the above measurements and 
upper limits in a halo-model independent fashion. We 
concentrate on light WIMPs with spin independent (SI) 
interactions, and extend the halo-independent method 
of Fox, Liu, and Weiner [TT], later extensively employed 
in [T2", by including the energy-dependent energy res- 
olution, efficiency, and form factors. In our form, the 
method can be used by any experiment to present their 
own results in a way that would allow for an immediate 
comparison between experiments in a halo-independent 



The differential recoil rate per unit detector mass, typ- 
ically in units of counts/kg/day/keV, for the scattering 
of WIMPs of mass m off nuclei of mass number A, atomic 



number Z , and mass z is 
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where E is the nucleus recoil energy, p is the local WIMP 
density, iia,z = 'rnmA,z I (fn -\- mA,z) is the WIMP- 
nucleus reduced mass, aA,z{E) is (a multiple of) the 
WIMP-nucleus differential cross-section daA.z /dE, and 
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is a velocity integral carrying the only dependence on the 
(time-dependent) distribution f{v,t) of WIMP velocities 
V relative to the detector. Here 



lmA,zE 
2f^A.z 



(3) 



is the minimum WIMP speed that can result in a recoil 
energy E in an elastic scattering with the A, Z nucleus. 
Due to the revolution of the Earth around the Sun, the 
77 function has an annual modulation generically well ap- 
proximated by the first terms of a harmonic series 

?7('^min, t) = 77o(Wmin) + 77l(''^min) COSUj{t - Iq) , (4) 

where lu — 27r/yr and to is the time of maximum signal. 

For spin-independent interactions (SI), the WIMP- 
nucleus cross-section can be written in terms of the ef- 
fective WIMP-neutron and WIMP-proton coupling con- 
stants /„ and fp as 

a'A'ziE)=<7p^[Z + iA-Z)ifJfp)r FlAE), (5) 

where Cp is the WIMP-proton cross-section and ^ (E) 
is a nuclear form factor, which we take to be a Helm 
form factor normalized to Fa,z{0) = 1. In most models 
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the couplings are isospin conserving, /„ = fp. Isospin- 
violating couplings /„ 7^ fp have been considered as a 
possibility to weaken the upper bounds obtained with 
heavier target elements, which being richer in neutrons 
than lighter elements, have their couplings to WIMPs 
suppressed for fn/fp — —0.7 [T3]. 

Fox, Liu, and Weiner [TT] observed that the factor 
'7(?^min) = o-p{p/m)r]{vrain) in Eq. ^ with SI interac- 
tions is common to all experiments, and compared direct 
detection experiments without any assumption about the 
dark halo of our galaxy by expressing the data in terms 
of Wmin and ?7(wmin)- This was done extensively in [T^, 
separately for ?7o(i'min) and 77i(t'inin)- Since the E-Vmhi 
relation depends explicitly on the WIMP mass m, this 
procedure can be carried out only by fixing m (except 
when m is much smaller than the masses of all nuclei 
involved, in which case the combination mvmin becomes 
independent of m). 

However most experiments do not measure the re- 
coil energy E directly, but rather a detected energy 
E' subject to measurement uncertainties and fluctua- 
tions. These are expressed in an energy response func- 
tion Ga,z{E, E') that incorporates the energy resolution 
aE{E') and the mean value (£") = E Qa.z{E), where 
Qa,z{E) is the quenching factor. In this context, recoil 
energies are often quoted in keVnr, while detected en- 
ergies are quoted in keVee (keV electron-equivalent) or 
directly in photoelectrons. Moreover, experiments have 
an overall counting efficiency or cut acceptance e(-E') that 
depends on A compound detector with mass fraction 
Ca z in nuclide A, Z has an expected event rate equal to 
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We observe that the factor fj{vmin) is common to all ex- 
periments also when the rates are expressed in terms of 
the detected energies E' as in Eq. This observation 
allows us to extend Fox et al.'s method to the more realis- 
tic case of flnite energy resolutions and i5'-dependent effi- 
ciencies, without restrictions on how rapidly these quan- 
tities change with energy. 

For this purpose, using dE — (4/i^ z/''^A,z)'>^min rf^minj 
we write the average of Eq. ([6| over a detected energy 
interval £'2] as 

R[E[,E'^]= j dv^ir,'R-[i'^^E'^]{Vm\n)fl{Vrain)- (7) 
^ 



Here we have defined the response function for SI WIMP 

-'yl.Z'^min/ 



interactions, with Ea,z = ^'^'^-^inl'^A.z, 



t:,si ( ^ _ 2 t;min Ca,z <^a,z{Ea,z) 
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When several energy bins are present, like when bin- 
ning the data in energy or computing the maximum gap 
upper limit, we label each energy interval with an index 



i and write Ri{t), TZ: 



SI I 



^), etc, for quantities belong- 



ing to the i-th energy interval. For example, binning the 
harmonic series in Eq. (H m energy gives 



Ri{t) ^ Rqi + Rii cos[u{t - to)] . 



(9) 



Our task is to gain knowledge on the functions r/o{vaun) 
and 77i(wmin) from measurements i?oi i Ai?oi and Ru ± 
ARii of Roi and Ru, respectively. This is possible 
when a range of detected energies [E[,E'2] corresponds 
to only one range of ?;inin values [wmin.i, 'yinin,2], for ex- 
ample when the measured rate is due to interactions 
with one nuclide only. In this case, [wmin.ii ''^min,2] is the 
Wniin interval where the response function 7?.^, ^/-^{vmin) 
is significantly different from zero. Ref. [H] approxi- 
mated this interval with fmin,i = Vmin{E[ — asiE'i)) and 
Wmin,2 = i^inin(£'2 -I" (Je{E'2)). When isotopcs of the same 
element are present, like for Xe or Ge, the "ymin inter- 
vals of the different isotopes almost completely overlap, 
and Wmin,i, Wmin,2 could bc the Cyi^z-weighted averages 
over the isotopes of the element. When there are nu- 
clides belonging to very different elements, like Ca and 
O in CRESST-II, a more complicated procedure should 
be followed (see below). 

Once the [E'i,E2] range has been mapped to a 
[wmin,i, w„iin,2] range, we can estimate the timin-weighted 
averages 
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In the case of binned data, these equations read 7701 — 
Rm/Af^ rhi = Rii/Af^, with errors = ARoJAf 
and Arju = ARu/Af'. 

Upper limits on binned data can be set by replacing 
R[E',E'] above with the measured upper limit. Upper 
limits on unbinned data can be set using the method of 
Fox et al. [TT]. The smallest decreasing function 77(w,„in) 
passing through a point {vg^ris) is the downward step 
function 77(wmin) — Vs for I'min < Vs and zero otherwise. 
The minimal event rate is correspondingly 



E', 



dE'GA^ziEA,z,E')e{E'). (8) 
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We use this equation in the maximum gap method [2] 
to bound the value of ris as a function of Vs for CDMS, 
XENONIO, XENONIOO, and SIMPLE unbinned data. 
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For compound detectors hke SIMPLE, Eq. ( 13 ) is equiv- 



alent but more transparent than the method in Appendix 
A.l of [H]. 

The data and detector properties we use are as fol- 
lows. (We acknowledge criticism of some experimental 
analyses [H], and try to be conservative.) 

CoGeNT. We use the list of events, quenching factor, 
efficiency, exposure times and cosmogenic background 
given in the 2011 CoGeNT data release [16]. We sep- 
arate the modulated and unmodulated parts with a chi- 
square fit after binning in energy and in 30-day time in- 
tervals (we fix the modulation phase to DAMA's best 
fit value of 152.5 days from January 1). We correct 
the unmodulated part by surface-event correction fac- 
tors C{E) = 1 — er^ /^c ^ which are similar to those in 
[m for Ec = 1.04 keVee ("CoGeNT high"), 0.92 keVee 
("CoGeNT med."), and 0.8 keVee ("CoGeNT low"). We 
leave it to the reader to subtract a possible constant back- 
ground contribution. 

CDMS. For the upper limit on the total event rate 
we use only the T1Z5 detector which gives the most 
stringent limits at low WIMP masses. The energy res- 
olution is [0.293^ + (0.056i;)2)]i/2^ and the range for 
the maximum gap method is 2 keV-20 keV. For the 
modulation amplitude we use the 95% upper bound of 
0.045 events/kg-day-keV for a modulation phase equal 
to DAMA's in the energy range 5 keV-11.9 keV [10]. 

DAM A. We read the modulation amplitudes from [T]. 
We consider scattering off Na only, since the I component 
is under threshold for low mass WIMPs and reasonable 
local Galactic escape velocity. We show results for two 
values of the Na quenching factor: 0.3 and 0.45 (the latter 
suggested in [H])- No channeling is included, as per [TO] . 

XENONIOO. The exposure is 48 kg x 100.9 days. 
We convert the energies of the three candidate events 
in Ref. [7] into SI values, and use the Poisson fiuctuation 
formula Eq. (15) in [2^ to compute the energy fiuctua- 
tions. We use the light efficiency function C^s in Fig. 1 
of [7] . We obtain the cut acceptance by multiplying two 
factors: the overall cut acceptance, which we set to a 
conservative value of 0.6 since it is unclear why in Fig. 2 
of [7] it would depend on the WIMP mass when expressed 
as a function of SI, and the S1/S2 discrimination accep- 
tance, taken from the just mentioned Fig. 2. We use 
a maximum gap method over the interval 4 < S'l < 30 
photoelectrons. 

XENONIO. We follow Ref. \^ and use only S2 without 
S1/S2 discrimination. The exposure is 1.2 kg xl2.5 days. 
We consider the 32 events within the 1.4 keV-10 keV 
acceptance box in the Phys. Rev. Lett, article (not the 
arxiv preprint, which had an S2 window cut). We take 
a conservative acceptance of 0.94. For the energy resolu- 
tion, we are more conservative than 0: we convert the 
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FIG. 1: Measurements and upper bounds on the unmod- 
ulated and modulated part of the velocity integral rj{vmin) 
as a function of iJmin. For this case of spin-independent 
isospin-symmetric couplings and WIMP mass of 9 GeV, the 
XENONIOO and CDMS bounds exclude all but the lowest 
energy CoGeNT and DAMA bins. 



quoted energies into number of electrons 7ie = EQy{E), 
with Qy{E) as in Eq. 1 of [B] with k = 0.11, and use the 
Poisson fluctuation formula in [21 . 

SIMPLE. We consider only Stage 2, with an exposure 
of 6.71 kg days and no observed candidate event. We 
take an efficiency rj^E) = 1 — exp{—r[{Edep/ Ethr) — 1]} 
with F — 4.2±0.3. With no events observed, the Poisson 
and maximum gap upper limits coincide. 

CRESS T-II. We take the histogram of events in Fig. 11 
of Ref. [3]. The acceptance is obtained by adding each 
module at its lower energy acceptance limit in their Ta- 
ble 1. The electromagnetic background is modeled as 
one e/7 event in the first energy bin of each module. 
The exposure is 730 kg days. We assume a maximum 
WIMP velocity in the Galaxy such that W recoils can 
be neglected. To take into account the Ca and O com- 
ponents, we follow the same philosophy as Method 2 in 
Appendix A. 2 of [12] . but, without having to assume a 
constant efficiency in each energy bin, we are able to 
cover the CRESST-II energy range without gaps with 
the following binning: three high-energy bins {i — 4, 5, 6) 
with scatterings off O only (assuming a maximum Wmin of 
- 750 km/s): [17, 20], [20, 23], and [23, 26] keV; and three 
corresponding low-energy bins (i = 1, 2, 3) with the same 
Winin range and scatterings off O and Ca: [11, 13], [13, 15], 
and [15, 17] keV. To avoid complications with the over- 
lap of the tails of the weight functions T^f ^(umin), we cut 
them outside the Vmin interval [?;,„in(_B^j), ?;niin(-E'2i)], i.e. 
we do not enlarge the Umin interval using the energy reso- 
lution. Having determined ?7oi = Roif-^oh for « = 4, 5, 6 
using O only in Aq^^, we estimate the Ca contribution 
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FIG. 2: Same as Fig. [l] but for spin-independent isospin- 
violating couplings /n//p ~ —0.7 and DAMA quenching fac- 
tor QNa = 0.45. The first CoGeNT energy bin and the 
first five DAMA energy bins are compatible with all existing 
bounds. The dashed lines show the expected »7o(fmin) (up- 
per line) and »)i(umin) (lower line) for a WIMP-proton cross 
section of 2 x 10~^^ cm^ and a Maxwellian halo with Solar 
System speed relative to the WIMPs of 220 km/s, Id veloc- 
ity dispersion 220/\/2 km/s, and Galactic escape speed 544 
km/s. 



expected 77o(wmin) (upper line) and 7?i(i'min) (lower line) 
for a WIMP-proton cross section of 2 x 10~^^ cm^ and 
a Maxwellian halo with Solar System speed relative to 
the WIMPs of 220 km/s. Id velocity dispersion 220/^2 
km/s, and WIMP Galactic escape speed 544 km/s. These 
theoretical lines, representing the canonical model of the 
dark halo with a low escape speed in accord with the 
most recent determinations |22j . pass tantalizingly close 
to the data points. Given our very limited knowledge of 
the velocity distribution of WIMPs in the Galactic halo, 
it is not inconceivable that a distribution function with 
a sharp cut-off at large velocities may account for the 
current data. It is therefore of the utmost interest that 
CDMS extend their modulation analysis to lower ener- 
gies, so as to confirm or exclude the spin-independent 
interpretation of the GoGeNT and DAMA annual mod- 
ulations over the full Vmin range. 

We thank Chris Savage for extensive discussions at the 
beginning of this work, and Peter Sorensen for providing 
us with the list of event energies in XENON 10 [S'. P.G. 
was supported in part by NSF grant PHY-1068111, and 
thanks KIAS for hospitality during part of his sabbatical 
year. G.G. was supported in part by DOE grant DE- 
FG03-91ER40662, Task C. 



to bins j = 1,2,3 as R^i^ = -^cLjVoj+s, where A^i^^ 
contains only Ca. Then to reduce the effect of the prop- 
agation of errors in subtracting the Ca contribution, we 
combine the three low-energy bins into one, obtaining for 

it ^ = EU(^oj - R'cij)/ EU 

Fig. [1] shows the measurements and upper bounds on 
fjo and ?7i for a 9 GeV WIMP with spin-independent 
isospin-symmetric couplings. Here we took QNa — 0.3. 
The XENONIOO and CDMS bounds exclude aU but the 
lowest energy GoGeNT and DAMA bins. Although a 
larger value of Qwa would shift the DAMA points upward 
and to the left by a factor (5Na/0.3, the tension between 
GoGeNT and DAMA on one side and XENONIOO and 
CDMS on the other is strong. Varying the WIMP mass 
from 6 to 12 GeV does not improve the situation. 

The tension is alleviated for isospin- violating couplings 
fn/ fp = —0.7, especially if the DAMA Na quenching fac- 
tor is taken as Qwa = 0-45 (see Fig. 2). In this case, the 
first GoGeNT energy bin and the first five DAMA en- 
ergy bins are compatible with all existing bounds. Since 
CDMS and CoGeNT both use Ge, the tension of the 
higher energy GoGeNT bins with the CDMS modulation 
constraint remains. If the unmodulated CoGeNT rate at 
high recoil energies is subtracted throughout the energy 
range, the relative modulation amplitude would have to 
be high, ~ 30%. 

As an intriguing possibility, the dashed lines show the 
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